#### Exercicio 3 ##

## Distncia entre cidades ##

distancias = c(NA,3949,3000,3927,3949,NA,1273,3188,3000,1273,NA,1827,3927,3188,1827,NA)
matdist=matrix(distancias, nrow=4, ncol=4)
cidades=list(c("Atenas", "Madri", "Paris", "Estocolmo"))
matdist=matrix(distancias, nrow=4, ncol=4, dimnames=cidades)
matdist=matrix(distancias, nrow=4, ncol=4, dimnames=list(c("Atenas", "Madri", "Paris", "Estocolmo"),c("Atenas", "Madri", "Paris", "Estocolmo")))
matdist[upper.tri(matdist, diag=TRUE)]=NA
resposta=as.dist(matdist, diag=FALSE, upper=FALSE)


## Criao de um data frame ##

cor=c("Claro","Claro","Claro","Escuro","Escuro","Escuro","Claro","Claro","Claro","Escuro","Escuro","Escuro","Claro","Claro","Claro","Escuro","Escuro","Escuro")
dieta=c("A", "A", "A","A", "A", "A","B", "B", "B","B", "B", "B","C","C","C","C","C","C")
var.peso=c(0.1,1.1,3.7,1.5,-0.1,2.0,5.7,-1.2,-1.5,0.6,-3.0,-0.3,3.0,-0.4,0.6,-0.2,0.3,1.5)
datf = data.frame(cor, dieta, var.peso)
tapply(datf$var.peso, INDEX=cor, FUN=mean)

## Criando uma Matriz ##

mat=matrix(data=rnorm(15, mean=10, sd=sqrt(3.6)), nrow=3, ncol=5)
mat
rownames(mat)=paste("L", 1:3)
colnames(mat)=paste("C", 1:5)
ColMean=c(apply(mat,1,mean),apply(mat,2,mean))
ColVar=c(apply(mat,1,var),apply(mat,2,var))
DataF=data.frame(ColMean,ColVar)

## Lendo e Salvando seus dados ###

St1_1_300=read.csv("St1 1 300 um_data.csv", header=TRUE,sep=",",)
names(St1_1_300)
table(St1_1_300$Class)
str(St1_1_300)
St1_1_300$Grandes = St1_1_300$Diameter..ESD.>=300 
table(St1_1_300$Class, St1_1_300$Grandes)
sp1=St1_1_300[St1_1_300$Class=="Carcass",]
Carcass.file=data.frame(sp1$Class, sp1$Diameter..ESD., sp1$Volume..ESD.)
help(save)
write.table(Carcass.file, quote=FALSE, sep="\t", file="Carcass.file.txt", row.names=FALSE)

## Classes de Objetos ##

iris
iris3
colnames(iris3)
names(iris)
class(iris)
class(iris3)
names(iris)=c("CompS","LargS","CompP", "LargP", "Espcies")
summary(iris)
### apply(as.matrix(iris, nrow=150,ncol=5), 2, FUN=mean) ### Minha dvida!
apply(iris3,2, FUN=mean)

## Acrescentando Dados de Sntese ##

arquivo=read.csv("esaligna.csv", header=TRUE,sep=",",)
head(arquivo)
summary(arquivo)
arquivo$SomaBio=arquivo$tronco + arquivo$folha
arquivo$areabasal=pi*(arquivo$talhao/2)^2
mean(arquivo$areabasal)
arquivo$arvores.grandes=arquivo$dap>10
arquivo2=arquivo[arquivo$arvores.grandes,]
write.table(arquivo2,sep="\t", file="arvores.grandes.txt", row.names=FALSE, quote=FALSE)



